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Abstract 

A complete procedure for identifying the area of convergence of 
blood drops originated from a single static source is presented. Both 
for bloodstains lying on an horizontal and on a vertical plane a com- 
plete study is developed, based on error analysis and on an oppor- 
tunely calculated planar joint probability density for the orientation 
of the horizontal projections of the tangent lines at the points of im- 
pact. The method is capable of generating a probabilistic map for the 
area of convergence, directly linking the angles of impact, and their 
uncertainties, to the projection on the ground of the point of origin 
of the drops. One of the objectives consists in providing a statistical 
definition of area of convergence, extending to this topic the mathe- 
matical accuracy of the calculation of the angle of impact in bloodstain 
pattern analysis (BPA). 

1 Introduction 

One of the most important contributions coming from the analysis of the 
bloodstains at a crime scene concerns the possibility of locating the area 
of origin of the projected blood. The drop colliding on an impact plane 
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generates a bloodstain that is typically elliptical in shape and from the ratio 
of the axes of this ellipse the angle of impact a can be estimated by mean 
of the well established formula a = arcsin ci/a, where d and a are the best 
estimations for the minor and major axis of the ellipse, respectively (see, 
G-g-, [D [2] and the pioneer works of Balthazard |3] and MacDonell [1]). The 
angle a refers to the angle between the plane of impact and the tangent to 
the trajectory at the point of impact. The collection of the tangent lines at 
the point of impact, for a certain group of drops, is commonly projected on 
the horizontal plane and, whenever the bloodstains are referred to the same 
production event, these projections will converge and cross in an area that 
overlaps the normal projection of the region of origin PQ The importance 
of this result, from the forensic point of view, cannot be overestimated. 

Looking the scene from top, the path followed by the flying drops develops 
in straight line^; as a direct consequence, the procedure for the identification 
of the area of convergence is essentially a geometrical matter and requires 
little, if any, recourse to physical concepts and calculations. It's just a matter 
of drawing lines and spotting intersections. 

This simple scheme is complicated by two main problems. First of all a 
mathematically satisfying definition of area of convergence is missing: the 
idea of area containing the intersections generated by lines ... [5] is not as 
precise as wished in this geometric context. The selection of the shape and 
of the dimensions of the area of convergence remains quite arbitrary and, 
apparently, it is not based on strict mathematical or geometrical considera- 
tions. 

The second point concerns the uncertainties on the measurements of the 
bloodstains, originating a propagated error in the definition of the tangents 
of impact [Bj. As a consequence of the analysis of the statistical errors, the 
intersections of the projections of the tangent lines at the point of impact 
will not result in single intersection points [TJ [8] , but more precisely in areas 
of probability of intersection. How should this be taken into account? 

The aim of this paper is to present a procedure that defines and identifies 
the area of convergence by mean of a statistical extension of the usual treat- 
ment of the uncertainties of the angles of impact of the bloodstains. The 
idea is to profit of the statistical errors in the estimation of the angles of 
impact to generate a planar probability map for the intersections. The areas 
of highest probability of intersection will be interpreted of highest 

^We neglect the unusual case where wind or other lateral forces act on the drops. 
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probability for the origin of the projection. Every area in the plane will be 
therefore uniquely associated to a number: the probability that the area con- 
tains the normal projection of the origin in the horizontal plane. In this way, 
the mathematical rigor that led to the estimation of the angles of impact 
will be transferred into the further steps of the analysis, within a coherent 
framework. The results of this procedure may also be the basis for subse- 
quent physical calculations of the trajectories of the drops and indicative for 
the quest for the point of origin. 

To achieve the goal, we will first recall in Sec. |2] the notions of tangent 
of impact, in the two different cases of stains deposited on horizontal and 
vertical surfaces. In that section we will also show how to combine the traces 
found in these different planes into a unique context of analysis for the deter- 
mination of the area of convergence. The core argument will be presented in 
Sec. [3] where, by mean of the definition of a suitable joint probability density 
for the orientation of the projection of the tangent of impact in the horizontal 
plane, a statistical treatment will be argued for calculating the probability 
that specific areas of the plane contain the normal projection in the ground 
of the point of origin. Practical applicative examples will be described in 
Sec. Hland discussed in Sec. |5l 



2 Angles of impact and uncertainties 

The analysis of the angles of impact, from a point of view that aims to the 
identification of the area of convergence, is quite different for stains lying on 
horizontal or on vertical surfaces. 

As a matter of fact, the elliptical bloodstains found on horizontal planes 
have the longest axis contained in the plane of the trajectory: the identifica- 
tion of the plane of the trajectory is immediate, as soon as one measures the 
angle of orientation of the stain, with respect to a certain reference. In this 



paper we will call this angle, with respect to the x axis, 7 (see Fig. 1(a)). 
We decide to define 7 in the range [0,27r[. The uncertainty related to the 
process of measure of 7 directly leads to the definition of a statistical error 
57, for this angle. 

If the bloodstain is lying on a vertical surface, instead, some little work is 
required. The tangent of impact is not in general contained in the same plane 



as one of the axes of the ellipse (see Fig. 1(b) ). To determine 7 (due to the 



presence of a vertical impact surface, 7 is here restricted in the range ]0,7r[) 
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(a) (b) 

Figure 1: The angle 7 between the x axis and the normal projection in the 
horizontal plane of the tangent of the ballistic trajectory of the blood drop 
at the point of collision. In (a) the case of a stain in the horizontal plane is 
shown: the measurement of 7 is immediate. In (b) the stain lies on a vertical 
plane: the angle 7 must be calculated from the dimensions of the ellipse 
and the angle (see Eq. [1]). We decided to define G ] — vr, vr], negative as 
depicted in (b). 



one must estimate the axes of the ellipse d and a and the angle G ] — vr, vr] 
between the vertical y axis and the major axis of the ellipse (our choice is 
done so to have the angle negative in Fig. 1(b)). Then one can easily 
calculate the angle 7 by employing the trigonometric formulaH: 



7 



arctan ■ 



sin d>\/ 1 
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7T + arctan ■ 
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. 7r/2 
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ii a ^ d and G ]0, 7r[ 
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(1) 



ii a = d 01 



0,7r 



The use of Eq. [T]for the determination of 7 can be the basis for the estimation 
of the error ^7 on 7. Following and extending the approach to the treatment 
of statistical error presented in [B], referring to the ISO 9000 standard [9], we 



^Eq. [l]can be easily demonstrated by just applying the definition of the trigonometric 
functions, and considering that the sine of the angle of impact equals the ratio of the axes 
of the ellipse: sin a = d/a. 
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may calculate the error ^7 as the square root of the variance of 7: 



57 



(2) 



where 5- is the uncertainty in the measure of " ■ " and where we have rea- 
sonably supposed that the covariance is zero, between the three variables. 
Passing through easy but cumbersome calculations, one finally obtains 



57 
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(3) 



for every |0| G ]0, 7r[ and d. For cj) = 11/2 one can easily check that Eq. [3] 
reduces to Eq. (21) of [6], in the case of zero covariance. Moreover, for — )■ 
then 57 — )■ (5 



(a^ — d'^)/d'^, and the dependence on the values of 6a and 5d 
is suppressed, as expected. The rapid increment of 57 for angles of impact 
around 7r/2, i.e. for a ~ is again an expected result, as in [0]. 

Despite its complexity, Eq. [3] has a direct applicability and can be easily 
embedded in any software for the bloodstain pattern analysi^. 

Concluding this section, we point out that, from this point of the analysis 
on, the statistical uncertainties for the angle 7 can be treated in the same 
way, both for stains on the horizontal and on the vertical planes. The only 
difference rests on how these uncertainties have been measured or calculated. 



3 A probability density map for the area of 
convergence 

The next step of the analysis consists in putting together the bloodstains, 
with their angles of impact and relative uncertainties, in order to detect a 
possible area of convergence. Let N be the number of drops (or, equivalently, 
stains) in analysis. After calculating the angles of impact, a collection of 

''The software AnTraGoS, realized by the author and employed by the forensic science 
service of the Italian Polizia di Stato, already embeds the algorithm for the calculation of 
^7 and of other angular uncertainties. 
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straight lines in the xz plane is at hand: these lines originate from the stains, 
and each one is oriented at a 7^ angle with respect to the x axis (with k = 
1, . . . , N). These lines are the best estimations for the normal projection of 
the trajectory on the horizontal plane. For sake of simplicity, from now on, we 
will refer to these lines as 7fc-lines. Each 7fc-line is determined with an angular 
uncertainty given by 6'~fk- According to the theory of random errors [TO] one 
can interpret the 7^ angles and their uncertainties in a statistical sense. The 
idea is then to probe the points in the xz plane in search of the areas of largest 
probability of intersection for the projections of the trajectories, given that 
their orientation can be associated to a certain statistical distribution. 



Referring to Fig. 2(a), we may imagine to pick a sample point P on 
the xz plane, which is supposed to be free of obstacles and obstructions, 
for simplicity. Connecting the point P with the planar projection of the 
point of impact a line is obtained, generating an angle 6k with the 7fc-line 
{9k E ] — 7r,7r]). Within the language of the stochastic processes applied to 
the error analysis [11], we can assert that the probability density function 
(PDF) associated to the passage of the projection of the trajectory of the 
drop k across the point P is described by a wrapped normal distribution, the 
analog of the Gaussian in circular statistics The form of this PDF is 



pdP) = pM = E -^P ( -^^^^ ) , (4) 



2^ 



with 0"^ = (57fc)^ playing the role of a variance in the zero-mean distribution 
of the variable 0^0. In the usual case where 57^ <^ vr, it is sufficient to 
consider the central term of the sum [i.e. the element with n = 0) and the 
following approximation holdil 



q2 



otherwise 



with Ak suitable normalization factor. 



^We denote here with capital letters the (stochastic) variables and with lower case 
letters the values assumed by these variables. 

^The validity of this approximation is evident if one realizes that in this case the 
contribution to the wrapped normal distribution coming from the tails of the original, 
unwrapped normal distribution is negligible. In the following, as an usual condition, we 
implicitly assume that 5^k ^ and adopt a PDF in the form of Eq. [H 
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(a) (b) 

Figure 2: The angle 9^ between the 7fc-hne and the hne connecting a generic 
point P of the xz plane and the projection in the xz plane of the point of 
impact. The only difference between the cases of deposits in horizontal (a) 
or in vertical (b) surfaces is that in the case (a) the point of impact and 
the projection of the point of impact itself coincide. In (a), for reasons of 
symmetry, the angle 6k is selected in the range ] — 7r,7r], and, if oriented as 
in the example figure, it is e.g. chosen to be negative. In (b), due to the 
presence of a vertical surface, additional limitations on the values of 9k are 
to be imposed: in this case 6k g] — 7a;, tt — 7/c[. 



These k symmetrical distributions are approximately Gaussian, except 
for the cutoff imposed at ±7r according to the limited range of the B^'s. 
The case illustrated in Fig. 2(b) can be treated in the same way but, due to 
the presence of an obstructing vertical surface of impact, the range of Ok in 
Eq. [5] is now restricted in the range ] — 7^, vr — 7fc[. Other restrictions may 
be here imposed to account for the presence of obstructions like furniture, 
walls, pillars, curtains, . . . 

Now the idea is to join the PDF's of all the stains: for every test point 
P of the xz plane, one may evaluate the 6'i, . . . , 6'^? angles and calculate the 
joint probability density function pjoint{P) = Pjoint{6i, ■ ■ ■ , On) relative to the 
event that the point P contains the horizontal projections of all the N stains. 
It is reasonable to assume the 0^ variables as independent and therefore, by 
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multiplying the PDF's, one finally gets 



TV ^ 

Pjoint{P) = A - W — exp 

k=l 




(6) 



where A is an overall normalization factor determined by integrating Pjoint{P) 
over all the possible and compatible points of the horizontal surface in the 
scene. 

The probability Prob (S) that a given area S contains the horizontal 
projection of the origin of the blood is then easily obtained by integrating 
this joint PDF over S: 



The integrals for the determination of the normalization factor A and 
in Eq. [7] can be easily solved, to an as good as wished approximation, by 
numerical algorithms. In practice, the estimation of the integrals can be 
limited to areas around the points of intersection of the 7fc-lines, where the 
convergence is expected to occur. This approximation is valid thanks to the 
rapid, exponential decay of the PDF's, far from the mean values. 

4 Practical use 

One of the most important advantages of the scheme presented in Sec. E] 
consists in the direct applicability of the method to an arbitrary large set of 
stains and in the easy realization of an automated procedure for the calcu- 
lation process, by mean of suitable algorithms. In this section we present an 
example of application, realized on a simple forensic software. 

First of all we notice that, even if the p— functions defined in Eq. Eland 
Eq. E] are continuous in their range of validity, a discretization procedure 
is necessary from a computational point of view. In the examples of this 
section, we have chosen, e.g., to sample the p^'s and of the p joint functions 
at every 5 cm, both in the x and in the z directions of the plane. The 
integrals are then calculated by supposing the p's constant in the so defined 

5 cm— squares. The sampling and the integrations can surely be made more 
accurate, but this level of approximation has proven to be sufficient in most 
real applicative situations. 
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Referring to Fig. [3] where different scales are used for convenience, the 
squares in purple cover areas where the probability of joint convergence ex- 
ceeds 10%, squares in orange are relative to probabilities in the range 1-10% 
and squares in green to probabilities of 0.5-1%. The 7fc-lines are in blue 
and the uncertainty in their orientation (±crfc) is in black. In the illustrative 
examples, the choice of the values of 7^ and 0"^ is arbitrary. 



The symmetric example in Fig. 3(a) refers to a simple case where the 



two stains have the same angular uncertainty. As expected, the calculated 
area of convergence is also symmetrical: the total probability associated to 
the colored area exceeds 95%, while the probability of convergence in the 
purple square, where the 7fc-lines intersect, is only about 11%. A slight non- 
symmetrical situation is instead depicted in Fig. 3(b), where, due to the 
smaller errors in the angles of impact, the resulting area of convergence is 
smaller (total probability >97%, probability in the purple square 15%). 



In Fig. 3(c) the effect of an increased distance of convergence is shown. 
The 5 cm— sampling is probably too fine at this scale and, even if the area 
of convergence is more regular in shape, it leaves out a lot of squares with 
probability less than 0.5%: as a result, the total probability of the colored 
area is only about 61%. 

An important example is presented in Fig. 3(d) here the three 7^ angles 
have more or less the same value. Even if the 7fc-lines converge in a rela- 
tively small area, the presence of an error in their determination results in 
a very large area of probability of convergence, oriented along the direction 
of these lines. The colored area contains only 48% of the total probability of 
convergence, indicating that a large error in the determination of the point 
of origin is to be expected. The example illustrates the problems related to 
the standard definition of area of convergence [S] and also highlights the im- 
portance of having differently oriented stains in the scene, to have a reduced 
overall error in the analysis (see, e.g., Sec. 3 of and [I3j). 

In any case, the validity of the assumptions of the method should be al- 



ways carefully verified, as shown in the problematic case of Fig. 3(e) , where 



the algorithm has calculated an area of convergence despite of the evident 
presence of stains belonging to different production events. The fact that the 
spotted area of convergence lies far from the effective area of intersection of 
the 7fc-lines can be indicative on an error in the assumption of single projec- 
tion. All the three angles for the stains to the left converge in a small area, 
while the stain to the right does supposedly belong to another projection: a 
clear indication that its study should be separated from the analysis of the 
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(e) (f) 

Figure 3: Some applicative examples of the probabilistic method for the 
identification of the area of convergence, as implemented in AnTraGoS, the 
software for the BPA of the forensic science office of the Italian Polizia di 
Stato. The color of the squares is related to the probability of joint conver- 
gence of the projected trajectories in the area: purple >10%, orange 1-10%, 
green 0.5-1%. See the text for a detailed description. 
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other three. 

The efficiency of the method and of its implementation are evident in the 
combined treatment of large numbers of drops (up to the order of 10^), as in 
the most common forensic cases. As soon as the angles are calculated with 
their errors, the computational time for the solution of the integral in Eq. [7] 
is practically zero and the area of convergence is immediately visualized. An 



example taken from a real case is depicted in Fig. 3(f) the purple area is 
representative of a probability of about 49%, while the total probability of 
the colored area amounts to 99%. 



5 Discussion 

In Sec. [2] we have presented a procedure for grouping the stains, deposited on 
horizontal and on vertical surfaces, into a single framework. The most im- 
portant result, derived from error analysis, is enclosed in Eq. |3l that provide 
a direct and easily embeddable formula for the estimation of the statistical 
uncertainty in the determination of the horizontal projection of the angle of 
impact. By extending the error analysis beyond the determination of the an- 
gle of impact, we could then construct in Sec. [3] a simple applicative scheme 
for a mathematical characterization of the area of convergence. 

The statistical approach can locate the areas of the horizontal plane where 
the highest levels of probability of convergence are concentrated. One suit- 
able definition of area of convergence may accordingly be stated: given a 
certain level of probability, it is, e.g., the smallest simply connected area 
where the local probabilities of convergence add up to that level. 

Different choices for the probability value spot areas that obviously may 
differ in shape or size, and a certain variability can also be ascribed to the 
arbitrariness on the discretization of the horizontal plane for the numerical 
estimation of the integral in Eq. [71 However, for all practical purposes, 
the accuracy of a discretization with 5 cm spacing is usually more than 
acceptable. 

This scheme, apart from being a natural extension of the error analysis 
in BPA, provides the expert with a tool for answering of the most impor- 
tant questions usually arising in a forensic context: what is the error in the 
determination of the area of convergence? The important issues and prob- 
lems connected to the standard definition of area of convergence, like the 
undetermined connection between the points of intersection and the areas 



11 



and the undervalued effect of the uncertainties, are easily overcome in this 



framework. As seen in the comment to Fig. 3(a), the small area where the 



7A;-lines intersect is effectively related to a small probability of convergence, 
and a wider area has to be considered to get higher values: an area that more 
or less covers the uncertainties, too. The application of the procedure to spe- 
cific situations like the one depicted in Fig. 3(d), where, due to the similar 
orientation of the angles of impact, the intersections of the 7fc-lines are not 
indicative for the convergence, also reveals the importance of the probabilis- 
tic approach: the area of convergence is effectively quite widespread, in that 
case. 

The stains with smaller errors are, as expected, the most relevant in 
the definition of the joint probability density of Eq. El great care should 
then be used in minimizing the errors in the previous steps of the BPA, 
i. e. in carefully selecting the stains, taking pictures and estimating the best 
fitting ellipses. Not surprisingly, large variances Sjk ultimately result in the 
uselessness of the drops in the overall analysis, especially in presence of many 
other stains with significantly lower variance. 

In this scheme, it's easy to include/exclude single stains in the analysis 
and comparatively check the effect of this inclusion/exclusion in the results: 
this recursive testing method can be of some help to identify the groups of 
stains with common origin. 

The results of repeated tests have also confirmed the validity of other 
previously published results in BPA: the errors in the determination of the 
area of convergence become larger when the drops are all in the same impact 
plane and when the 7-lines intersect with a very small angle ([71 IT3]). 

As seen in some of the examples, the area of convergence do not usually 
approximate a circular shape: this evidence is indicative of the necessity of a 
statistical approach that, like the one proposed, may reveal the asymmetries 
of the specific situation and may overcome the limits of single-valued error 
expressions like "the area of convergence is determined with an error of about 
20 cm" . 

The idealized situation sketched in this article can be improved by set- 
ting suitable contour conditions that account for the geometrical limitations 
imposed, e.g., by the presence of furniture, walls or other obstructions: ev- 
ery single drop might have its own limit values and the overall consistency 
of the hypotheses should be always carefully checked. Particular attention 
should be paid in the implementation and in the usage of algorithms that 
automatically generate a probability map for the area of convergence: the 
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results might not be precise enough to describe a specific situation, and in 
most manual check of the single PDF's is suggested. 

The further possible step of spotting the area/point of origin and its sta- 
tistical uncertainty appears as necessary as complex. In fact, even if the idea 
of pushing the error analysis into the field of the three dimensional recon- 
struction of the trajectory of the blood drop [H] seems more than reasonable 
at this point, some complications arise from the necessary introduction of 
other variables such as the air friction and the impact velocity of the drop: 
these physical quantities are generally estimated at crime scene with a rele- 
vant uncertainty [151 IIHl [HI [H] • Probably, for all practical purposes, a good 
estimation of the minimum and maximum for each linear dimension of the 
location of the point of origin can be obtained by repeatedly calculating the 
time-reverted equations of motion |T9l|8], selecting time by time different val- 
ues for the key physical quantities. However, our belief is that an extended 
error analysis should be proposed for this last stage of the BPA, too. In this 
way the general framework would certainly be more consistent and complete, 
as well as the forensic relevance of the results. 
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